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Abstract 

We calculate the next-to-leading order perturbative QCD corrections to the transverse momen- 
tum weighted single transverse spin asymmetry in Drell-Yan lepton pair production in hadronic 
collisions. We identify the splitting function relevant for the scale evolution of the twist-three 
quark-gluon correlation function. We comment on the consequences of our results for phenomenol- 
ogy. 
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I. INTRODUCTION 



Single transverse spin asymmetries (SSAs) in high energy hadronic reactions continue 
to attract much theoretical and experimental interest. They are defined as differences of 
cross sections when one of the initial hadrons' transverse spin is flipped, divided by the 
sum: An ~ {d(T{S±) — da{—S±))/ {da{S±) + da{—S±)). The theoretical description of SSAs 
ha. p.oven to be a challenge g =ince the leading coUinea. partonic cont.ibut.on to the 
asymmetries vanishes [2|]. Over the past few years, there have been a number of theoretical 
developments that have led to much progress in the exploration of the underlying physics 
for single spin asymmetry phenomena. These developments mainly follow two lines: the 
so-called transverse momentum dependent (TMD) approach , |lO| , which 

uses parton distributions and/or fragmentation functions that depend on partonic transverse 



momentum, and the twist-three quark-gluon correlation function approach 



15 



16 



13, 



Q, Q, Q, 



18|. More recently, it has been found in some cases that the two approaches are 



closely related and describe the same physics |l9l. l20[|. 

So far, however, phenomenological applications of the approaches have been limited to 
the "bare" parton model, that is, to the zeroth order of perturbation theory without any 
QCD corrections, as the latter were generally not available. This situation was remedied 



very recently when the leading-order (LO) kerne 



twist-three correlation functions were derived 



21 



s for the scale evolution of the relevant 
22l | . In this paper, we take a further step 



toward a more comprehensive QCD description of single-spin phenomena by calculating 
next-to-leading order (NLO) QCD corrections for a particular observable, the transverse 
momentum weighted SSA in Drell-Yan lepton pair production. 

As demonstrated by many examples, next-to-leading order perturbative QCD (pQCD) 
corrections are typically very important in hadronic processes. They often lead to signiflcant 
i^'-factors, and also allow estimates of the size of yet higher order corrections. Moreover, 
an NLO calculation for a particular physical process will provide a direct test of QCD 
factorization for the associated observable, complementing the general arguments for such 
factorizations [3, ^, [23]. One-loop pQCD corrections to the hard-scattering factors in the 
TMD factorization approach have been calculated for the observables associated with the 
so-called fc_|_-even TMD parton distributions [3]. For the related fc^-odd TMD observables, 
especially relevant for SSAs in various processes, there has not been any particular calcu- 
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lation so far. The same is true for the twist-three approach, where all calculations for the 



SSAs so far have been at leading order only 



15 



16| . Previous studies have shown that 



QCD corrections for higher-twist observables are much more complicated to obtain than 



for leading-twist ones 



24 



271]. This is related both to technical difficulties result- 



ing from more complex partonic states, and to possible mixing between higher-twist matrix 
elements 2J] . 



On the other hand, the recent developments, especially the consistency between the 
TMD approach and the twist-three approach found in |19|, have provided confidence in 
our understanding of the underlying theoretical description of single-spin phenomena. They 
naturally motivate a study of NLO corrections to SSA observables. The calculations and 



results of 



19| will be the starting point for our derivation of the NLO corrections to the SSA 



in the Drell-Yan process. In this process, a transversely polarized nucleon with momentum 
Pa scatters off an unpolarized nucleon (Pb) to produce a virtual photon with invariant mass 
Q and transverse momentum q±, which subsequently decays into a lepton pair. 



MPA,S^)p{PB)-^riQ',q±)+X 



+ X 



(1) 



where S± is the transverse polarization vector of the incident nucleon. The Drell-Yan pro- 
cess is bested suited as a first example for the calculation of NLO corrections to single-spin 
processes. It is among the simplest processes in hadronic scattering, and its single spin 
asymmetry is also kinematically simpler than those for other processes. For example, for 
scattering with single transverse polarization, the pair transverse momentum q± and the 
polarization vector S± are simply correlated as e^^S^q^ = l^^l |g_L| sin 0, where (p is the 
azimuthal angle of q± relative to that of S±. If q± is measured experimentally, the corre- 
sponding single spin asymmetry receives contributions from the so-called Sivers effect in the 
TMD approach, applicable when q± <^ Q, or from the twist-three Qiu-Sterman matrix ele- 
ments when q_i_ ~ Q. As we mentioned above, the two approaches coincide in the kinematic 



regime of overlap 



191]. In the following, we will make use of this fact. 



The transverse momentum of the virtual photon (or the lepton pair) generally depends 
on various transverse momenta in the process, namely those of the initial partons, and those 
generated by gluon radiation. For the cross section differential in transverse momentum, one 
has to be careful to classify the different contributions, as a TMD factorization only exists 



in the limit of small transverse momentum, q± ^ Q 



However, if we integrate over all 
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transverse momentum q±, the cross section will depend only on the longitudinal vcvovcveniMia 
fraction of the virtual photon, and a coUinear factorization approach will apply. For the 
single transverse-spin dependent cross section, we have to suitably weigh with transverse 
momentum in order to obtain a non- vanishing result, because the unintegrated cross section 



has linear dependence on q^. The weighted cross section is defined as 28l. |29[|: 



{q^Aa{S^)) ^ J d'qM ^^^^^^^^^ ' (2) 

where we have simplified the expression by omitting dependence on any other kinematic 
variables. Since the transverse momentum has been integrated out, the above weighted 



cross section can be properly formulated in the coUinear factorization approach [20], and 

can be factorized into parton distributions and/or twist-three correlation functions for the 

incident nucleons, and partonic hard-scattering functions. In case of the SSA, the quark- 

gluon correlation function in the polarized nucleon will be an important ingredient in the 

factorization formula. It will be part of the following factorization formula for the above 

g_i_-weighted cross section: 

d{q±A(r{S±)) _ f dxidx2dx\ 
dQ^ 

where cxo = 'iira'^^/SNcsQ'^, with s = {Pa + PbY, <i{,x') denotes the anti-quark distribution 
of the unpolarized nucleon, and Tp^q the Qiu-Sterman matrix element for quark q. We have 
restricted ourselves here to one quark flavor; extension to more flavors and to scattering off 
a quark from the unpolarized nucleon is trivial. In the following, we will drop the label q of 
Tp^q for simplicity. As indicated, Tp is a function of two separate light-cone variables, and 
thus the convolution over momentum fraction will include both, as we will see. Tp is defined 
as 



o-Q [— — ^Tp^qixi,X2)qix')n{xi,X2;x'), (3) 

J Xi X2 X' 



Tp{x,,X2) = f ^^i!Le^(-i^>-+(-2-.i)PlC-) 
J 47r 



X {PA,s\mmch''9Fa^icmc,v-)Hv-)\PA,s) , (4) 

where C is the proper gauge link to make the matrix element gauge invariant, and where 
the sums over color and spin indices are implicit. 

In the factorization formula Eq. (3), the hard-scattering function can be expanded as a 
series in the strong coupling constant, 

^ = ^{o) + ^^(i) + ... ^ (5) 

/TT 



where Ti^^^ is the leading order term, Ti^^^ the NLO one, and so forth. In the following, we 
will demonstrate that the above factorization formula is valid at NLO level. In particular, the 
collinear divergence can be factorized into the parton distribution and quark-gluon correla- 
tion function, whereas the hard coefficient function is free of any divergence. The real-gluon 
radiation diagrams have already been studied in 19|, and the results can be carried over to 
our present calculation with relatively little effort. We will compute the virtual corrections 
as well. It is important to check that the soft divergence in real-gluon radiation is canceled 
by that in the virtual diagrams, so that we are left with only collinear divergences, which 
can be absorbed into the parton distribution and/or the twist-three correlation function, 
where they give rise to the scale evolution of the distributions. 

The rest of the paper is organized as follows. In Sec. II, we will derive the leading order 
expression for the hard coefficient, and calculate the virtual correction at next-to- leading 
order. In Sec. Ill, we calculate the real-gluon radiation contributions, and combine them 
with the virtual corrections. We will show that the soft divergence is canceled in the sum, 
and that the remaining collinear divergence can be removed by collinear factorization. We 
conclude our paper in Sec. IV. 



II. BORN DIAGRAMS AND VIRTUAL CORRECTIONS 

At the leading order, the virtual photon is produced in the quark-antiquark annihilation 
subprocess. In order to obtain a non-vanishing weighted transverse-spin-dependent cross 
section (g_LAcr(5'^)), we have to include an initial state interaction as shown in Fig. 1, which 
provides the required phase [l^. We perform our calculations in covariant gauge. Let p' 
denote the momentum of the incident anti-quark, kqi that of the initial quark to the left 
of the cut, kq2 that on the right, and kg = kq2 — kgi the momentum of the polarized gluon 
attaching to the hard part. This attachment may take place on the left side of the cut, as 
shown in Fig. 1(a), or on the right side, as in 1(b). 

The polarized gluon is associated with a gauge potential A^, and one of the leading 
contributions comes from its component A'^. The gluon's momentum is dominated by 
XgP + kg^, where Xg is the longitudinal momentum fraction with respect to the polarized 
proton. The transverse momentum kg± flows through the perturbative diagram and returns 
to the polarized proton through the quark lines. The contribution to the single-transverse- 
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(a) 



FIG. 1: Leading order contribution to the weighted transverse spin- dependent cross section. 

spin asymmetry arises from terms linear in kg± which, when combined with yield d^A'^, 
a part of the gauge field strength tensor F^~^. In order to compute this contribution, we 
expand the partonic scattering amplitudes in terms of kg± up to the linear term. The 
weighted cross section can, in general, be written as 

d{q^Aa{S^)) _ e^'^Sj f d%, d%, k ■ O')) a(x')T (k k) (6) 

where Ta{kqi, kq2) represents the non-perturbative matrix element for the polarized nucleon 
with full momentum dependence on kqi and kq2. In the above weighted cross section, we 
have an explicit term in the integral, along with the hard partonic part H{kqi, kq2] Q^). 

To obtain the coUinear factorization formula Eq. (3), we have to perform a coUinear 
expansion of the hard partonic part. For the leading Born diagrams in Fig. 1, we find that 
q± is related to the transverse momenta of the two quark lines as: g_L = fcg2± for Fig- 1(a) 
and q± = kqi± for Fig. 1(b). Therefore, the contribution from Fig. 1(a) to the coUinear 
expansion of Eq. (6) will be 

^9 

where the propagator associated with the initial state interaction produces the pole at fc^ = 
0. The single spin asymmetry arises from the phase of this pole. Similarly, the contribution 
from Fig. 1(b) will be 



{qlH{kqi, kq2] g^)) |Fig.l{a) = T+N 'K^l. ' (7) 

\'^q2 Si/ 



{qlH{kq,, kq2, Q^)),,^,(,) = - _^k^^\^_^, Kr^ ■ (8) 



''q2 " kql) 

The total contribution is thus 



{qlH{kq,, kq2; Q'))Fig,,(,+,) - _ ^+ ^ _ - _ fc+ ) _ • 

(9) 



When integrated over the transverse and hght-cone-minus components of the two momenta 
kqi and kq2, the combined terms Ta{kgi, kq2) and kg± produce the matrix element Tr pjx, x). 
One then obtains the leading order contribution to the weighted cross section as 28|] 

^(^^(M = j xmmi - QVxx's) . (10) 

from which one can readily determine the leading order hard coefficient to be 

nixi,X2;x')=6il-—)5{l-z), (11) 

Xi 

where z = Q'^/s with s = Xix's. 

The above derivation shows that the Born kinematics greatly simplify the collinear ex- 
pansion for the hard partonic part, because of momentum conservation. We can utilize 
this feature in the calculations of the virtual corrections to the Born diagrams as well. At 
one- loop order, the virtual corrections contain two types of diagrams shown in Fig. 2. In 
the upper two diagrams (a, b) of Fig. 2 the polarized gluon attaches to the side of the cut 
opposite from the loop correction which is represented by a blob. On the other hand, in the 
lower two diagrams {c,d), the gluon attaches to the side that also has the loop. The loops 
are displayed in detail in Figs. 3 and 4. Fig. 3 is the usual vertex correction, plus self-energy 
diagrams. If the gluon is on the side of the loop, we have to attach the gluon to all possible 
places in the virtual diagrams, see Fig. 4. 

As before, the Born kinematics simplifies the collinear expansion in the calculations of 
these diagrams. For example, to obtain the contributions from Figs. 2(a,b), we can simply 
multiply the Born result in Eq. (9) by the known [s^ Drell-Yan virtual correction factor, 
giving 

{qlH{kq,,kq2;Q')) |F,g.2(a+b) = ) _ ^."^^^^ 

(12) 

where Cp = {N^ — l)/2Nc with Nc = 3 the number of colors. Here we have used dimensional 
regularization, with D = A — 2e space-time dimensions and /i the mass scale to be introduced 
in order to keep coupling constants dimensionless ^. In obtaining this result, it is essential 
that the one-loop virtual correction in the Drell-Yan process amounts to a simple multi- 
plicative factor to the vertex j'^. As for the Born diagram, the collinear expansion is trivial. 



2 3^2 

8 + TT^ 

e 



^ Note that we also follow Ref. t3fl] to absorb a factor (1 — e) in the normalization do, which universally 
appears in all matrix elements and hence does not affect the final results. 
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id) 



FIG. 2: One-loop virtual correction to the weighted cross section: the gluon attaches to the opposite 
side of the loop correction ( upper two diagrams ); the gluon attaches to the same side of the loop 
corrections (lower two diagrams) . 







FIG. 3: Corrections to the quark- antiquark-photon vertex, corresponding to the blob in the upper 
two diagrams of Fig. 



and the phase for the SSA comes from the initial state interaction, i.e., the denominator of 
Eq. (HID. 

The calculation of Figs. 2(c,d) is more cumbersome, but again the collinear expansion 
will receive contributions proportional to and kg2± just as in Eqs. ([7]), ([HD. After a 
lengthy calculation, we find that the diagrams in Fig. 4 lead to the following result for the 
initial state interaction contribution: 

//1-rr,.2\ 

Fig.2{c+d) - 



^9 



2 3 
l2 ~ 7 



(13) 
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FIG. 4: Corrections to the quark- antiquark-photon vertex with gluon attachment, corresponding to 
the blob in the lower two diagrams of Fig. \^ 



identical to that for Figs. 2(a,b). Substituting the resuhs in Eqs. ( 1T2|) and ( |T3l) into the 
coUinear expansion formula, and accounting for the leading-order phase space in D dimen- 
sions, we find the total virtual correction to the weighted single-spin cross section: 



I — e) J X x' 



3 
e 



(14) 



In the next section, we will calculate the real-gluon radiation contribution, and obtain the 
final result for the NLO correction. It is important to verify that the soft divergence in the 
above virtual corrections is canceled against that in the real diagrams. We will check this 
cancelation in the next section. 

III. REAL CORRECTIONS AND FINAL RESULTS 



The real-gluon emission contributions to the single-spin asymmetry for the Drell-Yan 



process have been computed in 



19|. In that paper, the focus was on the SSA at fixed q±, 



whereas in our present calculation we are considering the g^-weighted cross section, which 
involves integration over all q±. While we can still use much of the set-up of the calculations 
of 13] , we have to redo them in D = 4 — 2e dimensions. This is relatively straightforward. 
Another issue we need to address is the transverse-momentum flow in the diagrams. In 
the calculations performed in 1^, the virtual photon is an "observed" particle in the sense 



that its momentum is kept fixed in the colhnear expansion of the hard partonic scattering 
amplitudes H{kqi, kg2; Q^). Since we are integrating over all q± in the present calculation, one 
might think that the photon momentum could be a function of the transverse momenta of the 
initial partons and hence change the colhnear expansion. However, because of momentum 
conservation it turns out that it is sufficient to allow transverse-momentum flow only through 
the radiated gluon. 

Before going into the calculation, we note that we will only consider contributions asso- 
ciated with the twist-three quark-gluon correlation function. There are two types of such 
contributions: real-gluon emission in the subprocess (qg) + — ^ 7* + (?, and the quark-gluon 
Compton process (qg) + g ■y* + q. We do not consider in this study contributions entering 



32|. 



with a three-gluon twist-three correlation function |16l . I3l| . We also ignore contributions by 
the "axial" twist-three quark-gluon correlation function considered in Ref. 

The real-gluon radiation diagrams yield soft and colhnear divergences when integrated 
over the transverse momentum, which are regularized by dimensional regularization. Fol- 
lowing a standard procedure for the phase space integrals as in the case of the spin-averaged 



cross section 



we obtain the following expression for the real-gluon radiation contribution 



to the transverse momentum weighted spin-dependent cross section: 

d{q±Aa{S±)) as f dx dx' f Air^i^V 1 ^_2e f^, , 



dQ^ 2-K J X X' \ J r(l-c; JO 

d 

x—Tf{x,x) {Dqqq{x') + Dqgg{x')) 
+Tf{x,x) {Nl^q{x') + Nl^g{x')) 

+T^(x, X - X,) + } , (15) 

where v is related to the partonic center-of-mass scattering angle ^ by f = (1 + cos 9)/ 2. 
The above expression contains three contributions. The first two are the derivative and 
non- derivative terms from the soft-pole diagrams, respectively; the third is the contribution 
by hard-pole diagrams, which only have non-derivative pieces. In = 4 — 2e dimensions, 
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we obtain the following expressions for the corresponding partonic hard-scattering terms ^: 



D 



D 



1 -t 



u t 
- + - 

t u 



2Q^s 
iu 



2e 



N, 



qg 



is) 

m 



2(iV2 - 
1 1 



1 s 



:i-e)i A + ^ 

-t s 



2Q^u 
is 



+ 2e 



[Q 



qg 



2Nc -stu 
+e(Mp - us'^ + 2P + 3Ps - s^)] , 



2A/'c ■ ' s + u 



+ Tr 



s + u 



stu 

iQ^-if + Q^u^-e{s + u)is + i)s 
-sH 



(16) 
(17) 

(18) 

(19) 
(20) 

(21) 



,2(iV2-l) 

where s, t, and u are the partonic Mandelstam variables for the 2 — > 2 processes, which can 
be expressed in terms of and z as 



s = — , i = — — (1 — z){l — v), u = — — (1 — z)v . 

z z z 



(22) 



From the above expressions, we see that the integral over v will contain divergences when 
^ or ^ 1. The main task in this calculation is to separate these divergences and 
identify them as soft or collinear, so that they can be canceled appropriately. 

First, let us examine the derivative term from the soft-gluon pole contribution in the qq 
channel. After integrating over v, it becomes 



a, Avr/i^Y 1 1 fdxdx' jf 1\ „ , .^2, 

2^[-Q^) rxr^2iv^y-^^^"HwJ "^"'"^^ 

-zf 



d_ 

dx 



X-:^Tf{x, x) 



+ In. 



(23) 



where we have performed an integration by parts in order to simplify the 1/e term. The 
latter comes from a collinear divergence, which will be canceled by factorization into the 
evolved Tp function. There is no soft divergence in this term, which is expected because the 
virtual diagrams do not contribute to the derivative terms, as we saw earlier. 



^ Note that in contrast to 30| we average over the polarizations of the initial gluon in the qg subprocess by 
the factor 1/(2(1 — e)), as is customary in the MS scheme. 
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The soft-pole derivative terms in the qg channel give only finite contributions. It is easy 
to perform the phase space integration, and we obtain 



27r2(iV, 



11 



dx dx' 
X x' 



x^T,{x,x) 



g{x'){l + z') 



-Xl - z){A + Az'^ - 5z) 



(24) 



The non-derivative term from the soft-pole diagrams in the qq channel has both soft and 
coUinear divergences. After integrating over the phase space variable v, we get 



0-0 : 



27r V / r(l - e) 2N, 



L||!|:r,(.,x),-M{-|^(i-.)- 



1 z^ - 



TT 



+—(5(1 -z) + {z^ - 3z^ 



ln(l - z] 



Inz 

. — z 



{\-z) 



+ 



(25) 



The double-pole 1/e^ term represents a soft-coUinear divergence, which will eventually be 
canceled. The coUinear-divergent term oc 1/e will generate part of the splitting function 
for the evolution of the unpolarized quark distribution and/or the twist-three correlation 
function. 

Now we turn to the hard-pole contributions, which are only non-derivative. For these 
contributions, the two arguments in the twist-three quark-gluon correlation function are 
difi^erent and may depend on partonic kinematics. As a result, the v-integral is somewhat 
more involved, and for some parts the integral cannot be performed completely. First, we 
will separate these parts by introducing "plus" -distributions of the form 



'f+ 



dv 



9{v)-g{^) 



dv- 



= / dv 



1-v 



(26) 



in the integrand. The distributions arise from terms oc v ^,{l — v) ^ in the integrand which, 
when combined with the phase space factor {v{l — v))~^, give rise to identities of the form 



ln{v] 



+ 0{^) . 



(27) 



For the qq channel, the part that cannot be further integrated over v analytically then reads: 



Tp\x^x 



I 

x[(l-„(l-.))3 + .](^ + C,^-l-^) 



a. I dx dx' dv f 1 1 



l-v{l-z) 



q{x') 



(28) 



Both plus distributions are needed because the integrand is divergent at both i> — > and 
V ^ 1. The distributions guarantee that the integral over v is finite. Besides, as one can see. 
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there is no divergence in the hmit z ^ 1. The remaining part of the hard-pole contribution 
in the qq channel is then rather straightforward to obtain, and we obtain upon integration 
over V 



as 
2^ 



1 l + z 



ln(l - z) 



el 



z) + a + z) 



\nz 



+Tf(x, x) 



2N, 



■z + Cf 



;6{l-z) 



1 l+z' 



TX 



el 



+ {l + z' 



ln(l - z) 

1 - z 



+ 



In 2; 

1 - z 



-A 
(1- 



5(1 



(29) 



Again, we have both soft and collinear divergences. 

Similarly, for the hard-pole contribution in the qg channel, we have a regularized part 
that cannot be further integrated over v analytically: 

dx dx' dv 



X x' {1 — v)j^ 



f(l 



g{x 



X [{l-v{l- z)f + v'' {I- zfz] 



2(iV2 - 1; 



1 



(30) 



and a remaining singular part: 



as 
2^ 



-qT J (i_,)r(i_e) 



dx dx' 



X X' 



yTF{x,x)g{x') ly 2{N^-l) 



-z + Ti 



R 



l{z^ + (l-zf) + {z' + {l-zf)\n^-^ 



+ 1 



(31) 



There is no soft divergence. The collinear divergence will be canceled by factorization of the 
gluon sphtting contribution to the spin-averaged anti-quark distribution function. 

As we mentioned at the beginning, the soft divergence has to disappear after adding the 
contributions by the real-gluon radiation and virtual diagrams. Indeed this happens, as 
inspection of Eqs. f[T^ and fl2S]) . fl2^ shows. Specifically, the term oc 1/e^ from the soft- 
pole diagrams cancels that from the hard-pole diagrams associated with color-factor l/2Nc, 
and the remaining 1/e^ term from the hard-pole diagrams (associated with the color-factor 
Cp) cancels against that from the virtual diagrams. This is an important cross-check on the 
consistency of our calculations, and demonstrates the importance of the hard-pole diagrams. 

After cancelation of soft poles, the result will only contain collinear divergences. We find 
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for the remaining pole term 



as 
2^ 



diOc dtXj 
X x' 



Tp{x, x)q{x') 



2Cf 



l + z^\ ,^ 1 



+ {Tf{x,xz) - Tf{x, x))q{x'){CF + 
+TFix,x)gix')Tn{z^ + {1 - zY)} . 



1 A + z 



2N/1- z 



(32) 



The residue of this colhnear divergence contains the sphtting functions governing the evolu- 
tion of the anti-quark distribution in the unpolarized nucleon and the twist-three correlation 



function. For the former, we have 



q{x) = + g / ^ [<i{x')Cf + 9{x')Tn{z' + {1 - z) 



(33) 



where z = x/x' and q^^\x) denotes the "bare" leading order anti-quark distribution. Simi- 
larly, we obtain the collinear QCD correction to the Qiu-Sterman matrix element at equal 
momentum fractions: 



Tf(x,x) 



+ Cf + 





fdx' f 




~ V 




'l + z 


2Nc) 


_l-z 



-Tf{x', x'z) 



1 + z' 
1 - z 



— z 



Tf(x\x') 



(34) 



From this equation, we directly read off the scale evolution equation for the "diagonal" 
twist-three quark-gluon correlation function at Xi = X2 = x: 



d 



9 In /i 



2Tp(x,x;/i^) = ^^^^ j '^\^CfTf{x' ,x';iJ?) 



l+z' 



27r 

Nr 



(35) 



Tf{x', x'; fi'^)z — {Tf{x', x'; /i^) — Tf{x', x'z; /i^)) ~'~ 



2 

asin'^) f dx' 

271 



1 - z 



x' 



Tf{x',x';h ) . 



(36) 



This evolution equation could also have been obtained from the perturbative calculation 
of the quark Sivers function at large transverse momentum performed in 19|. However, 
we note that in that paper a boundary term (oc 6{1 — z)) in the derivative contribution 
was overlooked'^. After correcting for this term, the result of |l9| becomes consistent with 
that given above. Eq. (!35l) is also consistent with the results derived recently by different 



This term does not, however, affect the consistency of the twist-three and the TMD approaches estabUshed 



in 0. 
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methods [2l|, |22|. We note, however, that the evolution equations derived in [2l[] go beyond 
ours, as they also contain the contributions from additional operators such as three-gluon 
ones, which we are not considering here. 

A few comments on the evolution equation are in order. First, it is evident that the scale 
evolution of the "diagonal" (xi = X2 = x) function mixes with the function at Xi 7^ X2, 
implying that the equation as it stands is not closed. In other words, there will be a more 
general evolution equation for the full function Tf{xi,X2). This feature is quite general 



for higher-twist parton distributions and fragmentation functions [2j, |25!, |26|, |27|] . Second, 
there is no particular simplification of the evolution equation in the large- A'^c limit. This is 
different from what was discovered for the evolution equations for other twist-three quark 
distributions, such as hiix) and e{x), where the evolution equations are closed (diagonal) 



in that limit 



251]. However, we notice that the high-x part of the evolution equation, i.e.. 



the large z limit of the kernel in the integrand, is the same as that for the spin-averaged 
leading-twist quark distribution, because the term Cf{{1 + z^)/{\ — z))^ is the ordinary 
leading order quark splitting function. This property will have important phenomenological 



consequences 



21| for the behavior of the quark-gluon correlation function at high x and thus 



for SSAs in hadronic processes at forward angles. 

We note that the renormalization and evolution of general twist-three quark-gluon opera- 
tors has been extensively studied over the past two decades 



l26|, l27| . The above result 



24 



25 



for the evolution of the Qiu-Sterman matrix element should likely also be reproduced from 
the evolution equations discussed in some of these papers. However, we notice that the Qiu- 
Sterman matrix element corresponds to a very different projection of the general twist-three 
quark-gluon correlation function, and we do not expect that there will be a simple relation 
between the above evolution equation and those for other specific twist-three distributions 



such as Qf 



m 



24 



25 



x) . The comparison between the above results (or the ones of 2l|, |22]) and those 



is very important and will be addressed in the future. 



After MS subtraction of the coUinear divergences into the quark-gluon correlation function 
of the polarized nucleon and the anti-quark distribution of the unpolarized nucleon, we obtain 
the full NLO expression for the soft-gluon and hard-pole contributions^ to the transverse- 



^ We remind the reader that we do not consider contributions associated with soft-fermion poles or with a 
three-gluon twist-three correlation function. 
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momentum weighted single-spin dependent cross section in Drell-Yan lepton pair production 
in pp collisions: 

d{q^Aa{S^)) 



as f dx dx' ( , 2\ Ti 



27r J XX' 
1 . d 



In {CpVqq + Vqg^qg ® TF(a;, XZ] i/)) 
A* 



x^Tp(x, x; + In — ^ 



2Ar/ ax 



X {Cf{1 + z') + 



+ 2 



ln(l - ^) 



1-z 



In 2; 
l-z 



2z^ - 3;z2 _ 1 



1 



Tf{x,x;h ) + {-^ + Cf){1 + z)Tf{x,xz;h ) 



+Tf{x, x; h') (^{Cf + 2^) (1 - ^) + - 8)5(1 - z)^ 

+g{x'; n') [{Xg^TFix, x; //^))( ^^^, ^ )-(l - ^)(4 + 4^^ - 5^) 



+TF{x,x;fi')Tn{{z' + {l-zy)ln 



+Tf{x,x; /j!^ 



K 1 



27;^ + 48;^^ - 29z^) 



+ 2z{l-z) 



2(Ar2_i)6^ 

+ / ( ~ + T\ ^ ) ^^ T^ ^ ' /^^)^(^'; /^^) 

y 1-^V^^+ (1-^)+/ l-v{l-z) 



+ 



-Tf{x, X 



[1-v). 



l-v{l-z) 
1 



■,^i')g(x';^i') [(l-v(l-z)f + v\l-zrz] 



.2(iV2-l) ■ ^""l-vil-z). 
As expected, the logarithms containing the factorization scale enter with the splitting func- 
tions for the evolution of the twist-three quark-gluon correlation function and the twist-two 
anti-quark distribution. 

One important feature of this result is its behavior near "partonic threshold" , that is in 
the large- 2; limit of the integrand, corresponding to s ~ Q^, when the initial partons have 
"just enough" energy to produce the virtual photon. Setting the scale /i — Q, we have the 
following structure of the NLO correction in this case: 



(37) 



d{q^Aa{S^)) 
dQ^ 



Ois f dx dx . 2\-/ / 2\ 



(38) 



Here we have written out the "double-logarithmic" term which dominates near threshold in 
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the MS scheme. The eUipses denote terms that are subleading relative to this term. The 
structure of this expression is identical to that for the spin-averaged g^-integrated NLO cross 
section near threshold, 

da f dx dx' 



dQ 



2 



as f dx dx' , , I 9^ r ^ f\\i{l — z)\ 

^o:t^ / r1^x-ii^)q{x'-ii^) ^Cf[-: ^ 

2n J X x' I \ 1 — z J 



+ ... 



(39) 



This means that the soft gluon contribution is spin-independent. It contributes in the same 
way to the spin-averaged and single-spin-dependent cross sections, and will lead to the same 
soft-gluon threshold resummation effects to these cross sections, at least at the leading 
double logarithmic level. This observation is very similar to that made for the transverse 



momentum resummation in the Drell-Yan process 33|. 



IV. CONCLUSIONS 

In summary, we have derived the NLO perturbative-QCD correction to the transverse 
momentum weighted single spin asymmetry in Drell-Yan lepton pair production in hadronic 
collisions. In the calculation, we have shown that the collinear divergences can be absorbed 
into the NLO twist-three quark-gluon correlation function of the transversely polarized nu- 
cleon and the unpolarized quark distribution of the unpolarized nucleon. This procedure also 
determines the evolution equation for the "diagonal part" of the twist-three Qiu-Sterman 
matrix element at equal momentum fractions, Xi = X2- We have found this equation to be 



21 



221. 



consistent with the more complete one derived recently in Ref. 

Our calculations suggest that a general factorization formula (see Eq. (3)) exists for the 
transverse momentum weighted spin-dependent cross section in the Drell-Yan process, in 



extension of the general factorization arguments given in [2 

We have found that both the evolution kernel and the full NLO expression for the spin- 
dependent cross section become identical to their spin-averaged counterparts in the "thresh- 
old" limit s Q^, or z 1. This will likely have the phenomenological consequence 
that the single-spin asymmetry for the Drell-Yan process will be quite stable under NLO 
corrections, in particular when r = Q'^/s is large. 

It will be important to carry out further studies. We have mentioned already that it 
may be possible to derive the evolution of the twist-three correlation functions also from 
some of the results of [24^. Also, it will be important to derive NLO corrections also for 
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other processes. For example, extension to semi-inclusive deep inelastic scattering should 
be relatively straightforward to do. 
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